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We propose that back-to-back correlations in azimuthal angle of jets produced in collisions of 
unpolarized with transversely polarized proton beams could be used to determine Sivers functions. 
The corresponding single-spin asymmetry is not power-suppressed, but is subject to Sudakov sup- 
pression. We present estimates of the asymmetry (without and with Sudakov effects) for RHIC at 
jet transverse momenta of ~ 10 GeV and show that it may reach a few per cent or more and could 
provide access to the gluon Sivers function. 
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I. INTRODUCTION 



^SJ ' An important goal of ongoing experiments with polarized protons at BNL's Relativistic Heavy-Ion Collider 

J> I RHIC is to contribute to a better understanding of transverse-spin effects in QCD. Of particular interest are 

■ single-spin asymmetries An, obtained from scattering one transversely polarized proton off an unpolarized one. It 
' was found a long time ago in fixed-target experiments 0, Q that such p^p collisions can yield strongly asymmetric 

distributions of hadrons in the final state. The most famous examples are the sizable (0(10%)) An found in 
the process pp^ — > tt X, which express the fact that the produced pious have a preference to go to a particular 
, side of the plane spanned by the proton beam direction and the transverse spin direction. Recently, the STAR 

■ collaboration at RHIC has found that such large asymmetries persist even at collider energies It is fair 
to say that to date the asymmetries have defied a full understanding at the quark-gluon level in QCD. One 
reason for this is that in QCD An for a single- inclusive reaction is power-suppressed as l/p'^ in the hard scale 
given by the transverse momentum of the pion. This makes the formalism for describing the asymmetries 
rather complicated, compared to more standard leading-power observables in perturbative QCD. An attempt of 

^ ' an explanation for the observed An has been given within a formalism 4] that systematically treats the power- 
suppression of An in terms of higher-twist parton correlation functions. Alternatively, it has been proposed 
0, 0, 13, 13 that the dependence of parton distributions and fragmentations functions on a small "intrinsic" 
transverse momentum fc^ could be responsible for the asymmetries, through the interplay with the partonic 
; ^ ' elementary cross sections that are functions steeply falling with p-^. An is generated from the fc^-odd parts of 
' the partonic scatterings, which acquire an additional factor 1/p^, making the mechanism again effectively higher 
twist. Measurements of just An m pp^ — > tt X will not be sufficient to disentangle all these effects, and it has 
also been shown recently that for this power-suppressed observable the mechanisms are not all independent of 
one another "o]. 

There is, however, a class of observables for which the fc^-dependent distributions or fragmentation functions 
alone are relevant, and may actually lead to leading-power effects. These are observables directly sensitive to 
a small measured transverse momentum. In spin physics, the most well-known example in this class was given 
by Collins ^J. He proposed to consider the single-transverse spin asymmetry in semi- inclusive deeply-inelastic 
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scattering (SIDIS), ep ^ e'nX, where the pion is detected out of the scattering plane. This asymmetry may 
receive contributions from the fc"*" effects mentioned above: from the transverse momentum of the pion relative 
to its quark progenitor (the so-called Collins effect 0), or from the intrinsic k-^ of a parton in the initial proton 
(referred to as Sivers mechanism Q). Here we propose a new observable sensitive to the latter effect. 

More precisely, the Sivers effect is a correlation between the direction of the transverse spin of the proton 
and the transverse momentum direction of an unpolarized parton inside the proton |^]. The Sivers effect in the 
process pp^ — > tt X has first been analyzed in detail by Anselmino et al. 0, who extracted the Sivers functions 
for valence quarks from a fit to the data under the assumption that the asymmetry is solely due to this effect. 
Subsequently, the single spin asymmetry for the Drell-Yan process, which is another process that belongs to the 
class of "leading-power" observables mentioned above, was predicted [lfil |. 

Following in part the notation of Ref. [lo| . the number density of a parton f = u,u, . . . , g inside a proton with 
transverse polarization St and three-momentum P, is parameterized as 



fix, fe^, St) = fix, k^) + ^A^ fix, fc±) -^^-(-Pxfc^) 



(1) 



where is the quark's transverse momentum, with k-^ — \k^\. The function fix, k'^) is the unpolarized parton 
distribution, and A^/ denotes the Sivers function. As one can see, the correlation proposed by Sivers corresponds 
to a time- reversal odd triple product St ■ (-P x fc^). Since for stable initial hadrons no strong interaction phases 
are expected, until recently it was widely believed that the Sivers functions had to vanish identically. It was 
then discovered 0,0,0,0, however, that the time-reversal symmetry argument against the Sivers functions is 
invalidated by the presence of the Wilson lines in the operators defining the parton densities. These are required 
by gauge invariance and, as under time reversal future-pointing Wilson lines turn into past-pointing ones, the time 
reversal properties of the Sivers functions are non-trivial and permit them to be non- vanishing. It is intriguing that 
the possibility of a non- vanishing Sivers function emerges solely from the Wilson lines in QCD. Another aspect 
to the physics importance of the Sivers function is the fact that it arises as an interference of wave functions with 
angular momenta Jz = ±1/2 and hence contains information on parton orbital angular momentum |llL ll4L ITHI |. 

In this paper we will discuss a single-spin asymmetry in pp scattering that also belongs to the class of "leading- 
power" observables. The reaction we will consider is the inclusive production of jet pairs, pp^ jet^jeta A", for 
which the two jets are nearly back-to-back in azimuthal angle. This requirement makes the jet pairs sensitive 
to a small measured transverse momentum, and hence allows the single-spin asymmetry for the process to be of 
leading power. The asymmetry for this process should give direct access to the Sivers function. Actually, in 
contrast to SIDIS and to p^p — > nX, the observable we propose here has the feature that it is rather sensitive to 
the nonvalence contributions to the Sivers effect, in particular the gluon Sivers function. The latter has not been 
considered so far in any phenomenological asymmetry study. Its precise definition has been given by Mulders and 
Rodrigues (l^. They define a gluon correlation function (in the light-cone (^+ = 0) gauge) 



Mr^ix,kT)^ I ^^^^e^''<{P,S\F+\Q)F+^iO\P,S)y^o, 



(2^) 

which can be parameterized in terms of gluon distribution functions. In particular, one has 



(2) 



MTijix,kT)gr} = -xP 



Gix,kl) + 



M 



-Gt(x,4) 



(3) 



where G(a;,fcy) is the unpolarized transverse momentum dependent gluon distribution inside an unpolarized 
hadron. Eqs. (jSJ) and ||2J) can be viewed as an explicit definition of Eq. ^ in the gluon sector, i.e., Gt corresponds 
to gix, fc^), up to a k^ dependent normalization factor. 

As the above shows, the inclusion of transverse momentum dependent parton distributions is necessary when 
discussing Sivers effect asymmetries. This extension of the ordinary parton distributions that are functions of 
the lightcone momentum fractions only is not straightforward from a theoretical point of view. Factorization 
theorems involving fc"*" dependence are generally harder to derive. Also, Sudakov suppression effects become 
relevant. As is well-known from unpolarized hadron- hadron collisions, the average transverse momenta of pair 
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final states, naively associated with originating from intrinsic transverse momentum, are energy {\/s) dependent 
and can reach several GeV at collider energies [tR IT^ . Clearly, such large average transverse momenta are not 
to be attributed to intrinsic transverse momenta alone, but mostly to the transverse momentum broadening due 
to (soft) gluon emissions. In our study we will therefore include such Sudakov effects, albeit within a somewhat 
simplified treatment. We will use experimental data to obtain estimates for the average transverse momenta of 
initial and radiated partons in unpolarized hadron-hadron collisions, which we will subsequently use to obtain 
estimates for the Sivers effect asymmetry in the azimuthal angular distribution of jets with respect to opposite 
side jets. We will do this without and with inclusion of Sudakov factors; one has to keep in mind that in the 
first case the average "intrinsic" transverse momenta we will find will effectively contain a significant perturbative 
(Sudakov) component. It will be instructive to compare the results of the two analyses. 

Our study is also motivated by the favorable experimental situation at RHIC, where the STAR collaboration 
has recently presented data for a closely related back-to-back reaction. In the next section we will discuss 
back-to-back correlations in the unpolarized case and compare to the STAR data. In Sec. Ill we will then address 
the spin asymmetry in pp^ jet]^jet2^. Section IV presents our conclusions and a further discussion of some 
theoretical issues. 



II. JET CORRELATIONS IN UNPOLARIZED HADRON COLLISIONS 



We first consider the inclusive production of jet pairs in unpolarized proton-proton collisions, pp — > jetj^ jet2 X . 
We take each of the jets to have large transverse momentum. This implies the presence of short-distance phenom- 
ena, which may be separated from long-distance ones. More precisely, the cross section for this process factorizes 
into convolutions of parton distribution functions with partonic hard-scattering cross sections that may be eval- 
uated using QCD perturbation theory. To lowest order, the partonic subprocesses are the QCD two-parton to 
two-parton scatterings. If the cross section observable is defined in such a way that it is insensitive to transverse 
momenta of the initial partons or to particles radiated at small transverse momentum, the factorization is the 
standard "coUinear" one, and the convolutions are simply in terms of (light-cone) momentum fractions. The 
observable we are interested in is slightly more involved. We choose the two jets to be almost back-to-back when 
projected into the plane perpendicular to the direction of the beams, which is equivalent to the jets being sepa- 
rated by nearly A(j) = (jjj^ — (jjj-^ = tt in azimuth. Such a configuration directly corresponds to an observed small 
transverse momentum of the jet pair. In the case of two-by-two scattering of coUinear initial partons, the jets 
are exactly back-to-back in azimuth. Deviations from this may result from additional partons being radiated into 
the final state. If the jets have fairly large separation in azimuthal angle, the dominant contribution to the cross 
section will come from a single additional parton radiated into the final-state against which the two jets recoil. 
Closer to Acf) = tt, radiation is suppressed, and Sudakov effects become relevant^. Intrinsic transverse momenta 
of the initial partons may become important as well. 

In such nearly back-to-back situations, factorization is not necessarily lost; rather, a factorization theorem now 
needs to be formulated in terms of parton distributions depending on light-cone momentum fraction and transverse 
momentum. Factorization theorems of this type have been discussed for the simpler process e"*" e~ AB X, where 
A and B are two hadrons almost back-to-back [2(tI l2ll | , and for Drell-Yan type processes . For factorization 
to occur it is essential that transverse momenta of initial or radiated partons are linked to the observed small 
transverse momentum only kinematically, that is, by momentum conservation, but are neglected in the hard 
scattering. For pp jet]^jet2X, factorization at small measured transverse momentum of the pair has to our 
knowledge not yet been proven explicitly, but here we will assume that it falls in the class discussed in Ref. [2^ . 

In the following we first consider only intrinsic transverse momentum of the initial partons and neglect pertur- 
bative radiation of particles into the final state and Sudakov effects. The latter will generally be very relevant and 
we will include them afterwards to leading logarithm for the purpose of estimating them. A treatment beyond 



Hard three (or more) jet configurations may also contribute near A(p = n, if the additional hard parton happens to be almost in one 
plane with the two jets. Such situations are expected to be relatively rare because they are not associated with singular behavior 
of the perturbative cross section in the azimuthal back-to-back region. We ignore them in our analysis. 
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leading logarithm is fairly complicated for pp ^Jeti jetj X and not within the scope of this work. Our treatment 
follows the same strategy as applied in Refs. l2^. 

In the absence of additional partons being radiated into the final state, momentum conservation implies that 
the sum of the transverse momenta of the two jets is equal to the sum of the transverse momenta of the initial 
state partons. To be more explicit, we expect the cross section to be proportional to 

j cPkid'ki (ki + fc^ - Pj-^ - Pf^) hiki) h{ki) , (4) 

where and i-*^ are the transverse momenta of the jets, and the fi are the transverse momentum distributions 
of the initial partons. In general, the fi will also depend on the parton lightcone momentum fractions. The 

factor U may be thought of as a smeared-out 5'^ (-^ii ~'~ ) representing the standard transverse- momentum 
conservation for coUinear partons. 

Since the distributions fi are not known, we will assume Gaussian transverse momentum dependence for 
simplicity: 



Moreover, we will assume that the average transverse momentum squared is the same for all partons in the 
proton and independent of x, i.e. (fc^^) = (fc-^^) for z = 1,2. These simplifications may all be improved upon 
at a later stage, when there are data requiring a more sophisticated treatment, but here we will focus on the 
proof of principle rather than on making an accurate quantitative prediction. The assumption of Gaussians is for 
convenience and sufficient for our purpose. One obtains 

n {{ki^) + {ki^)) ' 

where = + Pj^ . One has 

(r^)2 ^Pl-^^ + P^^^ + 2P};Pl^ cos(A0) , (7) 
where A0 is the separation of the two jets in azimuth. Hence, one finds that 

U oze^'^"''°<^'^\ (8) 

where Kjj = Pj^Pj^l i^^"^) ■ This implies that U is peaked around Ai/i = vr as expected. Expansion for small 
5(j) = A()f) — TT shows a Gaussian behavior near the peak. 

Correlations in At/) for dijets have been studied in [2^. Measurements of A(/) distributions have also been 
performed both at the ISR j2^J, and in the fixed-target experiment E706 0,0], albeit not for two-jet correlations, 
but rather for pairs of nearly back-to-back leading hadrons, usually pions. Recently, the STAR collaboration 
at RHIC (BNL) has presented precise data on hadron A^ correlations from unpolarized proton-proton collisions 
at ^/S — 200 GeV (and from heavy ion collisions). Our study of spin effects in back-to-back reactions in the 
next section will be tailored to collisions at ^/S = 200 GeV at RHIC, so we will compare our approach to 
the STAR A(/) distribution data displayed for "same-sign" hadrons in —it/2 < Scj) < n/2 in Fig. In the 
STAR analysis, the first ("trigger") hadron was required to have 4 GeV < < 6 GeV, and the recoiling hadron 
2 GeV < P^ < P,^ . The pseudorapidities of both hadrons were within \rifii ^ \ < 0.7. As can be seen from Fig. ^ 
the peak at b(\) — At/) — tt = is clearly pronounced and appears consistent with a Gaussian behavior in b(^. 
One also notices that the distribution does not decrease to zero at large which could be indicative of the 
perturbative tail corresponding to hard-gluon emission. 

The above considerations for dijet production can be modified to apply to dihadron production. The first 
modification is that one needs to take into account fragmentation functions describing the probability with which 
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FIG. 1: Two-hadron azimuthal correlation in the back-to-back region. See Ref. 1§\ for details on the experimental definition 
of the correlation function. The data are from [l9l . and the solid curve is obtained using the distribution (|SJ and the widths 
in Eq. fitting the overall normalization and an additive constant. The dashed line shows the result of a fit to the data 
with inclusion of leading-logarithmic Sudakov effects. 

a final-state parton emerging from the hard scattering will yield the observed hadron. The hadron will take a 
light-cone momentum fraction z of the parton momentum p, i.e. in particular for a collinear fragmentation process: 



(9) 



for the transverse components, and hence r-^ — P^/zi + Pi/z2 in the formulas above. For our estimate we 
will use in the following only average values for the hadron transverse momenta over the experimental bins: 
i-Phi) ~ 4-5 GeV, {Pi) ~ 2.5 GeV. We have then determined the corresponding average Zi in the framework of a 
leading-order calculation of the unpolarized dihadron cross section in p p scattering at RHIC energy — 200 GeV 
and find (zi) « 0.45, (^2) ~ 0.25. Here we have used the CTEQ-5 |2g set of parton distribution functions and 
the fragmentation functions of 27], both at leading order. 

Additionally, there will be a transverse-momentum smearing also in the final state, that is, the observed hadron 
may be produced at some small transverse momentum relative to the parent parton ,28j , implying 



Pi = zivi + K , Pi = Z2pi + k 



2 I 



(10) 



with • fcj = 0. This will also have an influence on the Scf) distribution of the two produced hadrons. We 
estimate this effect by replacing (fej^^^ + {k^'^) in Eq. © by {k^'^) + {k^"^) + {k^'^)/{ziy + (fc^^)/(^2)^ where 
(ki'^) is the average transverse momentum broadening squared in fragmentation. We could improve this estimate 

by taking into account that for a final-state particle emitted at some angle only a certain projection of fc is 
relevant for the 5(j) distribution. This modifies the functional form of the distribution in 5(j). However, we found 
this effect to be rather insignificant for the 54) and the widths we consider below. 

Our next goal is to obtain an estimate of the average [ki ^) at RHIC energy = 200 GeV from a comparison 
to experimental data. For this, Ref. j2^ is particularly useful, where an analysis of a variety of data from 
fixed-target and collider experiments on transverse momenta of dimuon, diphoton, and dijet (or dihadron) pairs 
was performed. If one neglects radiative effects, such pair transverse momenta are directly related to intrinsic 
transverse momenta. The results of [29j | show that the pair transverse momenta increase with center-of-mass 
energy. Also, they are consistently larger for dihadron pairs than for diphoton or dimuon pairs. This may 
be understood from the presence of k^ smearing in fragmentation. Additionally, in contrast to diphotons or 
dimuons, dihadron cross sections are dominated by scatterings of initial gluons, which may have a somewhat 
larger k^ broadening. As mentioned above, we neglect this effect. From the results shown in j2^J we then 
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estimate for dihadrons at ^/S — 200 GeV 

+ {ki^) + {ki^)/{zi? + {ki^)/{z2? « 15 GeV^ , (11) 
and for non-fragmentation final states 

{ki^) + {k^^) w 9GeV2 . (12) 
From this, we estimate, assuming {k^'^) = (^2^^) and (k^^) — (fc^'^): 



(fc,^2)«2GeV, \J{k^^) K 0.5 GeY . (13) 

Our value for the width for the initial-state broadening is larger than that found from studies of single-particle 
inclusive cross sections at the lower fixed-target energies in flOj . The curve in Fig. ^ shows our result for the 
5(j) distribution based on the widths in Eq. (|13|l . We have fitted the overall normalization of the curve to the 
data, and we have also allowed an additive constant to the d(f> distribution in this fit, in order to account for the 
perturbative tail at larger \S(f>\. The resulting curve gives a fair description of the data, even though the data 
appear to prefer a somewhat larger width of the peak. 

The fairly large size of yjkf^ in Eq- l|13|l (as compared to typical hadronic mass scales) again indicates that 
there are significant perturbative effects that should be taken into account in a more thorough analysis. As we 
mentioned earlier, Sudakov effects, related to multi-soft-gluon emission, are expected to be particularly relevant. 
Near dcj) — 0, gluon radiation is kinematically inhibited, and the standard cancelations of infrared singularities 
between virtual and real diagrams lead to large logarithmic remainders in the partonic hard-scattering cross 
sections. For the S(j) distribution, these have the form ln^''~"^ {S(j))/S(j) in fcth order of perturbation theory, with 
1 < m < 2fc or, more generally, for the distribution they are of the form ln^''~™(s/|r^p)/|r-'^P (3Q> .sUl . 
where -s/l is the partonic center-of-mass energy. It is possible to resum these logarithmic contributions to all 
orders in as. For the leading logarithms (m = 1), this was achieved a long time ago Isol 133,^3. Recent progress 
in the resummation for — > AB X at next-to-leading logarithmic level was reported in 3^. Here we provide 
an estimate of the Sudakov effects by taking into account the tower of leading logarithms. 

Applying the derivation of |3^ to the process pp — > jetj^jetj X, the resummation of independent soft gluon 
emissions to leading double logarithmic order leads to the following distribution in r^: 

u'''"^-^ = TT J'^^'^^ \~- ^ ^^"ibM hib')Mb% (14) 

fJo « Jq Air I TT J 

where ao is the lowest-order cross section integrated over all r^, and C is the sum of the color charges for the 
external legs in the partonic hard scattering, i.e., for subprocesses involving quarks and antiquarks only one has 
C = ACf = 16/3, for processes with two quarks and two gluons C = 2{Ca + Cp) = 26/3, and for gg gg, 
C = AC A = 12. All these channels are relevant, because of the competition between the magnitude of the parton 
distributions and the magnitude of the Sudakov suppression. Finally, in Eq. 114|l 

Uh')= I d?ki^b-k^U{ki)^e-''(^^')'\ (15) 



where the last equality follows for our Gaussian k^ distributions of Eq. jSJ. In lowest order in a^, i.e. setting 
= in the Sudakov exponent, Eq. p4|l reduces to U of Eq. Q). 

We have used the Sudakov improved U^^^ of Eq. (fTljl in a fit of the value for (fc^^) to the STAR data. We 
find that inclusion of the leading logarithms leads to a markedly better agreement with the data, demonstrated 
by the dashed line in Fig. ^ and to a reasonably small value of (fc^ ^) « 0.9 GeV, which is closer to a typical 
hadronic mass scale and to the one obtained in |l(l |. 

We note that our treatment of the soft-gluon resummation for the 5(j) distribution somewhat differs from that 
developed in [33i. .35] , where the resummation was performed in terms of a one-dimensional integral transform. 
We found our Tpp^ach, which is based on a two-dimensional impact parameter h 0, to be numerically very 
similar to the one of [23 near (50 = 0, but to be better applicable out to larger 54>. A more complete study of 
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the Sudakov effects would start from the well-known |3l| full form of the Sudakov exponent, given in terms of 
an integral over gluon transverse momentum. This form reduces to ours in Eq. H14|l if the running of the strong 
coupling is neglected. For a running coupling, the Sudakov exponent becomes sensitive to the strong-coupling 
regime, since Us is also probed at scales near Aqcd- This will require the introduction of further non-perturbative 
contributions. It may also be important to determine next-to-leading logarithmic contributions to the exponent 
[33 |. The issue of matching to a fixed (next-to-leading) order calculation j^^] at larger Scf) will then become 
relevant. 

We finally note that, besides the peak around A(j) = tt, the two-hadron azimuthal correlation also shows a peak 
at A0 = 0, corresponding to the two hadrons being in the same jet. The width of this peak should be primarily 
related to the "fragmentation part" (k^ ^) / (zi)'^ -I- (fc^ ^) I (z-i)"^ of the width in Eq. One therefore expects the 

peak at A(/) = to be narrower than the one at A(/) = tt, which indeed it is for the STAR data. We have checked 
that a distribution of the form (jSj), with only fragmentation fc^ broadening as numerically given by Eqs. and 
l(T^. fits the = peak well. We note that a detailed description of this region will require using two-hadron 
fragmentation functions as studied in [^^i 1^ ■ In the following we will not consider further the same-side peak 
at A^ = 0. 



III. JET CORRELATIONS IN pp^ SCATTERING 

We will now study two-jet correlations near A^ = tt in the case that one of the two incoming protons is polarized 
transversely to its momentum. As we discussed in the introduction, our motivation is that such correlations may 
offer access to the Sivers function. The basic idea is very simple. As follows from Eq. the Sivers function 
represents a correlation of the form St ■ (P x fc"*" ) between the transverse proton polarization vector, its momentum, 
and the transverse momentum of the parton relative to the proton direction. In other words, if there is a Sivers- 
type correlation then there will be a preference for partons to have a component of intrinsic transverse momentum 
to one side, perpendicular to both St and P. Suppose now for simplicity that one observes a jet in the direction 
of the proton polarization vector, as shown in Fig. [5] A "left-right" imbalance in of the parton will then affect 
the A(j) distribution of jets nearly opposite to the first jet and give the cross section an asymmetric piece around 
A(j) = TT. The spin asymmetry 

^ J (16) 

will extract this piece. 

The denominator of this asymmetry will be, up to normalization, the function U discussed in Sec. II. We now 
define the y, z directions as given by the polarization and momentum, respectively, of the polarized proton. The 
numerator may then be found by considering Eq. Q with an additional factor k^ ^ in the integrand, corresponding 
to the Sivers correlation in the polarized proton: 



V 



[ <fiki<fki 5^ (ki + Pf^ Pf \ h{ki) hikt) , (17) 



where we have introduced the bar on fi to indicate that /i is related to the transverse momentum dependent 
part of a Sivers function, with Gaussian width {k^"^). The latter will in general be different from that in the 
unpolarized distribution (in fact, smaller to satisfy a positivity bound, see Ref. We have normalized k^^ 

by Vi'^i^) instead of by \k^\ as Eq. ^ suggests. This follows the analysis of Ref. ^3 (d- Eq- below) and 
takes care of the fact that for k^ = the Sivers effect should vanish. 
In this way we obtain 



n{{.i^) + {ki^))' 
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FIG. 2: Asymmetric jet correlation. The proton beams run perpendicular to the drawing, 
or in terms of the jet transverse momenta and azimuthal angles (measured w.r.t. St), 



(\Pf^\ sin cf^ J, + \Pf J sin (I) J, 



) ^ 



As an example, for the case of 4>j-^ — 0, corresponding to our specific example displayed in Fig. |3 this yields 



V cx sin(A0)e 



-K-p cos(A0) 



where = 2PI^PI:J {{ni^) + {ki^)). 

Setting {k'^'^) — r(fc^^), the resulting spin asymmetry is 



An = 



(20) 



(21) 



where for the denominator of this asymmetry we have taken the function U given in Eq. (jSJ with the assumption 
{k:t^) = (k^^) for z = 1,2. 

There are, however, several different partonic channels contributing to jet production in pp scattering. For 
each of these, there is a combination of parton distributions, with dependence on light-cone momentum fraction 
and intrinsic transverse momentum. Thus in general, there will be a weighted sum of fc^ dependent functions in 
the numerator and the denominator of the asymmetry, that is. 



(22) 



Here, the A''^/i(xi, k^) are the Sivers functions as introduced in Eq. Q. The convolutions ® are over light-cone 
momentum fractions only. We note that since the Sivers functions correspond to distributions of unpolarized 
partons, the hard-scattering cross sections o'/j/j are the usual unpolarized ones in both the numerator and the 
denominator of the asymmetry. The pseudorapidities of the jets are denoted by rjj-^ and rjj^ . The hard scattering 
functions depend only on large scales, that is, on Pj[ and Pj^. Therefore, any dependence on fej^j is neglected in 
the (T/i/j; in other words, one considers the first term in a collinear expansion. 

One can see that a simple result such as Eq. (|21|l will only emerge if the x and k^ dependences in all functions 
factorize from each other, if all distributions in the numerator and, separately, in the denominator depend on 
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fc^ in the same way, and if the x-dependence for each Sivers function is identical to that of the corresponding 
unpolarized one. 

We win now give simple estimates for the possible size of the spin asymmetry A-^. For this purpose we will 
need a model for the dependence of the parton distributions on the light-cone momentum fraction x as well as 
on transverse momentum. We will assume that the x and dependences may indeed be separated for each 
function. For the unpolarized densities we write: 

/(.,fc^)=/(.)-^e-('=^)V<'=^^), (23) 

where the f{x) are the usual unpolarized light-cone distributions for flavors f = u,u, . . . , g, for which we again 
use the CTEQ-5 leading order set For the moment, we will neglect Sudakov effects, so we will use the 

value -y/ {k-^ ^) = 2 GeV we found in the previous section for the initial-state broadening without rcsummation of 
Sudakov logarithms. 

The parameterizations for the Sivers function we will use are taken from jl^ , where they were determined from 
comparisons to data of Ref. fl^ on inclusive single spin asymmetries for p^p — > ttX: 

A^/(x,fc^) = 2AA,(x)/(a:)— ^v/2^(T^-^e-('=")^/<«"^> . (24) 

We have chosen this parameterization and the accompanying fit values of Ref. ^3 since these are the only ones 
available so far. The question of universality (i.e. process independence) of the Sivers function we leave as an 
unresolved issue. For a further discussion of this point see Sec. IV. 

As before, r — {k-^ ^) / {k-^ ^) , and we use r = 0.7, in accordance with the fit of Ref. j^. This value allows for a 
good fit to the E704 data QJI and also to A polarization data, as discussed in Ref. 0| . The A// (x) are x-dependent 
normalizations, defined in [l^ as 



(at- + bfY^f+^'f') 



af^fbf f 

In Ref. |lCl| |. only the valence u and d Sivers functions were taken into account, since the data of fl| are in the 
forward region of the polarized proton, corresponding to large momentum fractions in its parton distribution 
functions. For u and rf, the parameters extracted from comparison to the data read ,.1 Q] : 

Nu - 0.5, au = 2.0, 6„ = 0.3 , 

Nd - -1.0, aa = 1.5, = 0.2 . (26) 

For the sea quarks we assume relations identical to 125|) . with Nu{x) = J\fu{x), N^g{x) = Afd{x) for simplicity. 
The details of these choices are not crucial. 

The size of the asymmetry is, however, very sensitive to the gluon Sivers function. The reason for this is that 
in the kinematic regime we will investigate here, ^/s — 200 GeV, Pj^ ^ 10 GeV, contributions from gluon-gluon 
and quark- gluon scattering are most important, because of their large partonic cross sections and because the 
parton momentum fractions become as low as xi,2 ''^ 0.05. There is so far no experimental information on the 
gluon Sivers function. It may be possible to obtain estimates within models of nucleon structure. This could be 
an interesting topic for future studies, but is beyond the scope of the present paper. To illustrate the dependence 
of the asymmetry An on the gluon Sivers function, wc will simply present results for four distinct cases: 

(i) Afg{x) = {Mu{x)+Md{x))/2, 

(ii) Ng{x) = , 

(iii) Mg{x) ^Nd{x) , 



(iv) Ng{x) = Nd{x) , but with = 2.5 GeV for gluons. 
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FIG. 3: Predictions for the spin asymmetry An for back-to-back dijet production at RHIC (-^s = 200 GeV), for various 
different models (i)-(iv) for the gluon Sivers function (see text). The solid line marked as "(iii)-l-Sud." shows the impact 
of leading logarithmic Sudakov effects on the asymmetry for model (iii) . 

These choices represent a very small (scenario ii), an average (with respect to the functions for quarks, scenario i), 
and two somewhat larger (scenarios iii and iv) gluon Sivers functions. Other choices are clearly possible, including 
yet larger functions. Choice (iv) is motivated by the notion that intrinsic smearing could be larger for gluons. 
We keep r = (k^^)/ {k^'^) = 0.7 in all cases, although this choice does not result from the fit of Ref. where 
a gluon Sivers function was not included. We have found that a change in r for the gluons of 10 % leads to a 
change of about 6 % in the asymmetry at its peak. 

The resulting asymmetries, as functions of 54), are shown in Fig. |21 at RHIC's ^/s = 200 GcV. For simplicity, 
we have chosen <f>j-^ = 0. We have taken into account jets with pseudorapiditics jvy^j 2I — 1 (^^ suitable for the 
STAR detector) and 8 GeV < i^^J;, < 12 GeV. The strong sensitivity to the gluon Sivers function is evident from 
Fig. IHl One can see that sizable asymmetries are by all means possible. In fact, the asymmetry may easily be 
even much larger (> 10%) if the gluon Sivers function is favorably close to the unpolarized gluon density. We 
expect asymmetries of ^ 1% to be easily measurable at RHIC. A typical value for the statistical uncertainty of 
such measurements may be estimated as [4ll | 

« , (27) 

where P is the beam polarization, C the integrated luminosity, and a the unpolarized dijet cross section integrated 
over the kinematical bin we are interested in. Using P — 0.5, a moderate C — 10/pb, and estimating a — 6-10^ pb, 
we find SAj^ « 2 • 10^*. It will of course be important to understand systematical uncertainties at a similar level. 



It is straightforward to determine the angle where the asymmetry has its maximum. If we choose (pj-^ — and 
define (5(/> = (pj.^ — tt, then we find: 

cos(^0.„..) « 1 - , (28) 

which as expected is a function of the observed jet transverse momenta, and of (k^^) and (fc^^). For our 
parameters given above, this yields 50max ~ 0.48 (for scenarios i, ii and iii). The value of the asymmetry at this 
S(f>ma.yi dcpcnds of coursc on the magnitude of the Sivers effect functions. 

As in the previous section we will also estimate the effect of Sudakov factors by including soft gluon emissions 
at the leading double logarithmic level. Their effect on the denominator U of the asymmetry has been described 
in Eq. H14|l . A difference is now that for jet pairs (unlike the case of inclusive hadron pairs) only the initial-state 
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broadening plays a role j42l |. This means that we will now have C = 8/3 for processes with only initial quarks 
and/or antiquaries, C = 13/3 for qg scattering, and C = 6 for a gg initial state. For the numerator of the 
asymmetry one has to replace in Eq. H14|) 

h{h^) - j Skie^b-^^^^hiki) = l6-e-^^<«^^)/4 ^\bi,{h^) . (29) 

This leads to 

V = J{n^^) / — 6Ji(6r^) exp C h{b')f2{b') . (30) 

The resulting spin asymmetry is also shown in Fig. |3 for the case (iii) above, i.e., J^g{x) = Md{x). For this curve 
we have now for consistency used the smaller value {k^ 2) w 0.9 GeV that we extracted from our Sudakov 
analysis of the STAR data in the previous section. This leads to a shift of the peak of the distribution closer 
to 5<j) = 0. More importantly, it is evident that the inclusion of the Sudakov factors leads to a considerable 
suppression of An- Of course, this does not rule out a sizable spin asymmetry in principle. As we mentioned 
earlier on, depending on the normalization of the gluon Sivers function, we could have a much larger than given 
by models (i)-(iv). In this context, we would also like to remark that the spin asymmetry for the single-inclusive 
reaction pp^ it X at moderately high p-^ , from which in 10] the valence quark Sivers functions were extracted, 
is also sensitive to a small transverse momentum, and hence susceptible to Sudakov effects. The analysis of Ref. 
flOj did not include Sudakov factors, which means that any effects of Sudakov suppression have effectively been 
absorbed into the Sivers function itself. In that sense the curve in Fig. 13 may include Sudakov suppression more 
than once. A Sudakov improved analysis of the asymmetry in pp^ ^ -k X would therefore be desirable. 

To gain statistics in experiment, one will not just select events with (pj-^ « and vary (pj^] rather one will want 
to integrate over bins in (pj-^ . Here one has to take care not to wash out the asymmetry by simply integrating over 
all (f>j-^ . The asymmetry will in general have the following dependence on and S(j): 

v4n((?!'ji , (5(/)) cx \Pj'_^\sin(pj^ — {Pj-J {sin pj^ cos S(j) + cos (pj^ sin S(j)) A{cos{S4})). (31) 

One possibility is to select "jet 1" in the hemisphere of the "spin-up" direction, and to weight the asymmetry 
with cos(0jj ) over this hemisphere: 

A]S!{S(p) = / d(j}j^ cos{4) j^) A]S!{(pjj^,S(j)) (X sin S(j)A{cos{S(p)). (32) 



One has to keep in mind that this weighted asymmetry has a complicated dependence on the transverse momenta 
of the two jets, which differ per event. In general, it will be rather involved to extract the normalization and the 
width of the Sivers functions from such a weighted asymmetry, but the above is one of the ways to obtain a Sivers 
asymmetry that is dependent on only the angle Scj). More generally, one could project out the full ^ azimuthal 
angular dependence, preferably for a fine binning in the transverse momenta of the jets. 

If one were to consider leading hadrons instead of the jets, one would run into the problem that there could be 
additional effects generating a single spin asymmetry at leading power. In the fragmentation process the Collins 
effect could contribute 0, which is a correlation between the transverse spin of a fragmenting quark and the 
transverse momentum direction of the outgoing hadron relative to that quark. We do not, however, think that 
this mechanism will be very important here. First of all, the fragmenting quark would need to have inherited 
its spin from the transverse spin of the proton. This means that the transversity densities of the proton would 
be involved, and that the partonic cross sections would depend on transverse spin in the initial and final states. 
These cross sections are much smaller than the unpolarized ones we used in our study above In addition, a 
major difference between the Collins and Sivers effects is that the gluon Sivers function is allowed to be nonzero, 
whereas a gluonic Collins functions is forbidden by helicity conservation. We are therefore confident that studies 
ofpp^ ^ hi h2 X (with hi, /12 two hadrons almost back-to-back in azimuth) at RHIC, for example in the PHENIX 
experiment (where it would complement the Drell-Yan single spin asymmetry measurements) , would also be useful 
for learning about the Sivers functions. 
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IV. CONCLUDING REMARKS 

We have proposed an observable that could provide access to the Sivers effect: a single transverse-spin asym- 
metry in the distribution in relative azimuthal angle « tt of jets in a dijet pair. Unlike the more customary 
single-spin asymmetries for single- inclusive final states in pp^ scattering, this observable is not power-suppressed 
in a large energy scale. It also has the advantage that it will be directly (and only) sensitive to the Sivers func- 
tions, in contrast to An for the process pp^ —> tt X for which several different competing mechanisms could be at 
work. Using experimental information on the average transverse momentum of initial and radiated partons from 
dimuon, diphoton, dijet and dihadron production in hadron-hadron collisions and using results from the study of 
Ref. m on the Sivers effect in pp^ ^ tt X, we have presented estimates for this new observable. These indicate 
that the asymmetry could well be at the few percent level, which should make it experimentally accessible at 
RHIC. 

Our further analysis revealed, however, that Sudakov effects lead to a significant suppression of the asymmetry. 
We stress that this does not necessarily mean that the asymmetry must be small. It turns out that the unknown 
gluon Sivers function mainly drives the size of the asymmetry. We know of no theoretical reason why this 
distribution function should be small. In any case, any sign in experiment of the asymmetry we propose will be 
definitive evidence for the Sivers effect. We also point out that the distribution in azimuthal angle between the 
jets is only one example of a variety of similar observables in pp scattering. Other closely related examples, which 
deserve further attention and may be equally suited for experimental studies, are the total transverse momentum 
of the jet pair, the distribution in "Pout" (the momentum of one jet out of the plane spanned by the beam axis 
and the other jet's transverse momentum) [s^, and the distribution in "P^-balance" z — —Pf^ ■ E3- 
Any of these may be obtained from different projections in the two-dimensional transverse-momentum plane and 
may be predicted using our formulas above. 

We close with a few comments on some theoretical issues, which we hope also provide directions for future 
work. As we mentioned earlier on, for observables that have a hard scale but additionally involve an observed 
small transverse momentum, factorization theorems are rather hard to establish. For our back-to-back dijet 
distribution, the issue of whether or not factorization occurs still remains to be investigated. For the case of 
nontrivial polarization effects, such as the Sivers or Collins effects, this is a particularly relevant issue since, unlike 
in the unpolarized case, the effect itself already relies on the presence of non-perturbative "intrinsic" transverse 
momentum: the parton distribution functions need to have an intrinsically nonperturbative dependence on the 
transverse momentum, arising as jM or k-^ /{k^"^). 

In case factorization does apply, another related complication is the apparent non-universality of the Sivers 
functions. When it was recognized that the presence of Wilson lines allows the Sivers functions to be non- 
vanishing 1^ 0, IT^ IT^j l , also the remarkable result followed that the Sivers function in SIDIS differs by a sign 
from the one that enters in the Drell-Yan process. This process dependence is a unique prediction of QCD. 
It is entirely calculable, but has not been studied yet for other processes, such as pp^ — > ttX or the reaction 
pp^ — > 2 jets X we have considered here. The color gauge invariant factorization is expected to make definitive 
statements here. Jet reactions may generally be easier to analyze theoretically than reactions with observed 
hadrons in the final state, since the latter involve also fragmentation functions that inevitably complicate the 
issue of gauge invariance further 9, ,i44 | . A novel aspect in all this will also be the color gauge invariant definition 
of the gluon Sivers function, which has so far not been obtained (the same applies to any transverse momentum 
dependent gluon distribution and fragmentation functions). It will result in the proper gauge invariant version of 
Eq. ^ . For our present study, the process dependence implies an uncertainty. We have refrained from making any 
ad hoc choices and simply used as a starting point the valence quark Sivers functions obtained from pp^ tt X 

m- 

It is evident that it will be very important to deal with these issues. A proof of a factorization theorem for 
the process pp — > 2 jets X at small pair transverse momentum would be highly desirable. That said, even if 
factorization will be shown not to occur, the observable we have proposed is obviously still a quantity of interest. 
It will presumably then give us insight into novel aspects of QCD dynamics. 
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